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Abstract NMR quantum information processing stud- 
ies rely on the reconstruction of the density matrix rep- 
resenting the so-called pseudo-pure states (PPS). An 
initially pure part of a PPS state undergoes unitary 
and non-unitary (relaxation) transformations during a 
computation process, causing a "loss of purity" until the 
equilibrium is reached. Besides, upon relaxation, the 
nuclear polarization varies in time, a fact which must 
be taken into account when comparing density matrices 
at different instants. Attempting to use time- fixed nor- 
malization procedures when relaxation is present, leads 
to various anomalies on matrices populations. On this 
paper we propose a method which takes into account 
the time-dependence of the normalization factor. From 
a generic form for the deviation density matrix an ex- 
pression for the relaxing initial pure state is deduced. 
The method is exemplified with an experiment of relax- 
ation of the concurrence of a pseudo-entangled state, 
which exhibits the phenomenon of sudden death, and 
the relaxation of the Wigner function of a pseudo-cat 
state. 



A. Gavini-Viana ■ R. S. Sarthour • I. S. Oliveira 

Centro Brasileiro de Pesquisas Ffsicas, Rua Dr. Xavier Sigaud 

150, Rio de Janeiro 22290-180, RJ, Brazil 

E-mail: gavini@cbpf.br 

A. M. Souza 

Institute for Quantum Computing and Department of Physics 
and Astronomy, University of Waterloo, Waterloo, Ontario, N2L 
3G1, Canada. 

D. O. Soares-Pinto ■ E. R. deAzevedo • T. J. Bonagamba 
Instituto de Fi'sica de Sao Carlos, Univcrsidade de Sao Paulo, 
P.O. Box 369, Sao Carlos 13560-970, SP, Brazil 

J. Teles 

Instituto de Ciencias Exatas e da Terra, Campus do Medio 
Araguaia, Universidade Federal de Mato Grosso, Rodovia MT- 
100, Pontal do Araguaia 78698-000, MT, Brazil 



Keywords Pseudo-pure states • NMR Relaxation 
Entanglement dynamics ■ Dccoherence 



1 Introduction 

The dynamics of open quantum systems has been used 
for a long time to study physical ideas behind quan- 
tum measurements and decoherence pQ. Special atten- 
tion has been devoted to those aspects of open quantum 
evolution which are relevant to understand the transi- 
tion from quantum to classical regimes. Recently, the 
need for the maintenance of quantum state coherence, 
for quantum computation and communication tasks [2] , 
has increased the interest in dynamical properties of 
open quantum systems. In general, relaxation occurs as 
a result of the coupling of the system with the environ- 
ment, which leads to loss of coherence, being charac- 
terized mainly by energy relaxation and phase random- 
ization rates [3] . Energy dissipation carries away energy 
quanta on the scale of the Larmor frequency by enabling 
the spins to flip amongst their energy levels so as to es- 
tablish the equilibrium populations differences. Phase 
randomization destroys the coherences due to inhomo- 
geneities and fluctuations in the environment. Further- 
more, for nuclear spin systems with / > 1/2, additional 
quadrupolar interactions also leads to relaxation when 
modulated by the environment. Recently, theoretical as 
well as experimental studies about the nature of the 
spin-environment coupling have been reported [H[51[5]. 

On the other hand, Nuclear Magnetic Resonance 
(NMR) has been successfully used as an experimen- 
tal method for many Quantum Information Processing 
(QIP) implementations 7 . Generally, the NMR density 
matrix is written in the form: 
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Pe = (1 - + epi 

where p/ is the (normalized) identity matrix, and p\ is 
an arbitrary density matrix [5] . e is a parameter varying 
between and 1, which is basically the ratio between 
the magnetic and thermal energies. For liquid samples 
at room temperature, e ps 10 -5 , what makes p c a highly 
mixed state, close to the maximum entropy one. 

In most NMR QIP implementations, the experiment 
is carried out in an ensemble of spins in such a highly 
mixed state. Upon a set of unitary and non-unitary 
operations, the density matrix p\ can be transformed 
into a pure state, \i/j)(ip\. In this situation, p e is called 
a pseudo-pure state |10j . 

Pseudo-pure states allow exploring ensemble quan- 
tum computing using NMR technique because, upon 
unitary transformations, which form the basis of quan- 
tum logic gates and algorithms, only the pure deviation 
matrix evolves. The identification of the pure part of the 
pseudo-pure density matrix can be made from quantum 
state tomography technique (QST). 

However, PPS states are non-equilibrium states, and 
they relax back to a Boltzmann distribution in a NMR 
characteristic time Ti, the longitudinal (spin-lattice) 
relaxation time. Besides, NMR coherences vanish in a 
time T2, the transverse relaxation time. Regarding re- 
laxation studies of pseudo-pure states in the context of 
NMR quantum information processing, one can point 
out the studies reported on references [TTUT5lfT51[T5] . In 
those studies, relaxation behavior of various pseudo- 
pure states have been investigated in different NMR 
systems. In reference it was found that spin-lattice 
relaxation followed multi-exponential model that in- 
cluded mixed magnetic dipolar and electric quadrupo- 
lar interactions. In Ref. [12] it was found that whereas 
cross-correlations accelerate the relaxation of certain 
pseudo-pure states, they delay that of others. On the 
other hand, a robust method for quantum process to- 
mography provided a set of Lindblad operators that op- 
timally fit the density matrices measured at a sequence 
of time points [13] . 

In spite of the various studies cited above, to the 
best of our knowledge, there is no report showing how 
the initial pure part evolves under relaxation processes. 
Yet, it is just this part of the density matrix that mat- 
ters for NMR QIP analysis. The reason is that there is 
no simple normalization procedure to keep the devia- 
tion part of p e a true density matrix along the relax- 
ation. Therefore, it would be desirable to have a way to 
follow the evolution of an initial PPS state which un- 
dergoes a relaxation process \ip)(ip\ — > pi, yet keeping 
the properties of a density matrix for the deviation, all 



the way through. This would allow the study of other 
important quantum information observables under re- 
laxation. In this paper we propose such a method. To il- 
lustrate it, we analyze experimental data for the concur- 
rence of initially pseudo-entangled states and the dis- 
crete Wigner functions of a relaxing quadrupolar spin 
3/2. 

The paper is organized as follows: In Sec[2]we show 
that there is no trivial normalization procedure for the 
deviation density matrix, and that fixed normalization 
schemes produces invalid density matrices. In Secf3]we 
develop a normalization procedure that produces valid 
relaxing density matrices. Finally, in Secf4] we present 
some application examples of the method. 

2 Time-Independent Normalizing procedure for 
the deviation density matrix 

The density matrix formalism is an appropriate descrip- 
tion for both, pure and mixed quantum states [2]. Here, 
we consider a general density matrix for a spin ensemble 
with the form: 

p = p r + a (1) 

where pi = 1/N represents a maximally mixed state 
and a is called deviation density matrix [9]. This type 
of density matrix appears, for instance, from spin en- 
semble coupled to a heat reservoir at a fixed tempera- 
ture and whose thermal energy is much bigger than the 
magnetic energy of nuclear moments in the ensemble. 
This is called the high-temperature approximation. 

For some systems with zero trace observables, such 
as the case of NMR, only a can be experimentally ob- 
served. This means that quantum information quanti- 
ties, such as the concurrence, entropy or Wigner func- 
tions, among others, cannot be directly obtained from 
experimental data; for that a true density matrix is nec- 
essary 

We consider now a quantum process described by a 
trace preserving quantum operation S. Applying such 
operator to the density matrix Eq.([T]) yields 

S(p)=S(p I )+S(a) (2) 
= p I + [S(p I )-p I + S(a)], (3) 

which has the same form as ([T]) , if we make the identi- 
fication: 

a^[S( Pl )- Pl }+S(a), (4) 
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This shows that the process iS is not generally linear in 
<t. The linearity of S with respect to a is only obtained 
if S is a unital process (5(1) = 1), for instance, unitary 
transformations, which is not the case for most relax- 
ation phenomena |14j . Trace-preserving unital processes 
are the quantum analogs of the classical doubly stochas- 
tic processes, which can only increase the von-Neumann 
entropy of a state (defined to be — Trjplogp}). Other- 
wise, non-unital processes provide the only means to 
reduce entropy and prepare initial states [Hj. Usually 
in quantum algorithms applications, the identity term 
of the Eq.JT]) can be neglected because it is invariant 
under unital processes. 

In order to illustrate a relaxation process S acting 
on the density matrix, the deviation density matrix ob- 
tained from Redfield theory for a purely quadrupolar 
relaxation of an ensemble of spin 7 = 3/2 will be con- 
sidered (the deviation density matrix elements are on 
Appendix |A"|) . We consider as initial state a pseudo-pure 
state 



Pe = {l-e) Pl + e\^){ijj\, 



(5) 



where \tp) is a pure state. It could be, for instance, the 
pure state |11). Prom §5$ we have: 



|ll)(ll| = -[p e -(l-e)pz] 
e 

1 

= -a + Pl . 

€ 



(0) 
(7) 



Hence, normalizing the deviation density matrix by e 
and adding to pi the pure state can be obtained. But, 
since a changes with evolving time due to quadrupo- 
lar relaxation process, the loss of purity will make the 
pure state |11) to change to a mixed state. In Fig. [1] 
(top) is shown the relaxation of the population (diag- 
onal matrix elements) of the |11) state. This was ob- 
tained considering fixed normalization, like Eq.([7| for 
er varying in time. Since populations assume negative 
values at times above 10 ms, one no longer has a valid 
density matrix. Therefore, maintaining the normaliza- 
tion of the Eq.0, leads to matrices exhibiting negative 
populations at some instants of time. On another hand, 
if we try adjusting the final population values, the equi- 
librium state p eqi we can use the following normaliza- 
tion 



(bottom) that they have non-negative values for the 
whole interval of time, but at t = we have incorrect 
populations values for the pseudo-pure state. Thus we 
also discard this method of normalization. 

From these two different ways to normalize initial 
and final deviation matrices, we can conclude that the 
normalization must be time-dependent. 



3 Time-Dependent Normalization Method 

First we will consider initial and final conditions for 
p(t). Re- writing (5) in a slightly different way, we have 
the pseudo-pure state, 



p(t = 0)=p I + e(\ip}(iP\-p I ). 



(9) 



In words, the pseudo-pure state is composed by a max- 
imum statistical mixture, pi, and a pure contribution 

Under relaxation, the system reaches the equilib- 
rium state. For this state one can consider the Zeeman 
interaction as the most important one. In the high- 
temperature limit: 



p(t — > oo) 



Pi 



H(3uj 



Tr{e^ 1 ') r± N 
pi + 2Ie (p z - Pl ) , 



(10) 



where / is the angular momentum of the nuclear spin 
and e = hf3uj/2, and p z = (1/NI)I Z + Pl . 

From this, we can write for the respective initial and 
final (equilibrium) deviation matrices: ctq = e (|^)(^)| — pi) 
and (Too = 2Ie(p z — pj). The initial deviation matrix 
has "polarization" e, whereas the equilibrium has "po- 
larization" 21 e. The relaxation phenomenon connects 
these two matrices in time. Therefore, we can write a 
general form for a deviation density matrix: 



a (p a - pi) , 



(11) 



where a is the polarization and p a the respective den- 
sity matrix. From a and pi, p a can be determined for 
any deviation density matrix a. 

Adding ao to both sides of the Eq. (jlip and isolating 
p a , we have: 



Peq = -0- + P I (8) 

such that the p eq has trace equal 1 and non-negative 
populations. In this case (electric quadrupole relaxation 
to 3/2 nuclear spin) e' = 3e. It can be seen in Fig. Q] 



p a = -\ 1 p)(^\ + -[a-a ]+(l--)pi. (12) 
a a \ a/ 

This is the final expression for the density matrix for a 
given a and a. For example, at t = 0, a = e and a = ao, 
and we are left with a pure state. Upon relaxation the 
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other terms come into action. Notice that the first and 
last terms on the right side can be interpreted as a depo- 
larization channel [2] , whereas the middle term depends 
on the relaxation model, for instance, the Redfield the- 
ory, a is a time-adjustable parameter, which varies in 
the interval (0, 2Je]. It is worth mentioning that there is 
a possible connection between a and the effective num- 
ber of nuclei per state [10] . 

Eq. (JT3J) can be used to obtain the properly normal- 
ized density matrix p ai at any instant of time, from 
the measured deviation density matrices. A numerical 
example is shown on Fig. [2J for the same parameters 
used to calculate the data on Fig. [TJ At any instant of 
time the populations add 1 and never become negative. 
Next we will analyze experimental results for pseudo- 
entangled states and Wigner functions. 



4 Experimental results 

23 Na NMR experiments were performed using a mag- 
netic field oriented lyotropic liquid-crystal system (Sodium 
Dodecyl Sulfate, SDS), using a 9.4 T - VARIAN INOVA 
spectrometer. The sample composition was 21.3% of 
SDS, 3.6% of decanol, and 75.1% of deuterium oxide. 
The quadrupolar coupling was found to be (16700 ± 
70) Hz at 24 °C. In the experiments that characterized 
the relaxation of all elements of the density matrix, an 
initial pseudo-pure state was prepared using the SMP 
technique [16j . In this case, all the deviation density ma- 
trix elements were measured using the QST method via 
coherence selection Q7] . The basic experimental scheme 
consisted of: i) a state preparation period performed 
with SMP; ii) a variable evolution period r where re- 
laxation takes place; and iii) a hard RF pulse with the 
correct phase cycling and duration to execute QST via 
coherence selection, as described in detail in Ref. [18] . 



4.1 Computational basis states 

In Fig. [3] we show some a-polarization curves for differ- 
ent initial states, namely the computational basis states 
1 00) , 1 01) , 1 1 0) and |11), obtained from experimental re- 
laxation results [18] . 

The density matrix for a computational basis state 
is diagonal and the "polarization" for each of these 
states is different. This happens because the relaxation 
acts differently on each diagonal elements of the density 
matrix. Since each state population relax at a different 
rate [11], the normalization will also be different. 



4.2 Pseudo- entangled state 

Entanglement is a quantum correlation that reveals non- 
local properties of a given system. In order to compute 
such correlation or, in other words, to measure the de- 
gree of entanglement within the system, it is possible to 
use the concept of concurrence [20] . For a given density 
matrix p, the concurrence can be written as: 

£(/)) = max {o, y% - yfia - v% - v 7 ^} , (13) 

where the A^s are the eigenvalues, in decreasing order, 
of the matrix R = p{(j y ®Uy)p* (uy®a v ). It is important 
to note that for density matrices like 

"a OOiu' 



_w* d_ 

named X— state, the concurrence has a analytical solu- 
tion given by: 

C(p) = 2max{0,C / ;C // } , (15) 

where C 1 = \z\ - Va~d e C 11 = \w\ - VbZ [2T][22]. 

For a NMR system, described by Eq.©, it was 
shown that even being \ip) a pure entangled state, i.e., 
a Bell basis state, the total density matrix is always 
separable [8j. It happens because the polarization is 
very small (e 10 -5 , as said before), meaning that 
the density matrix is highly mixed. However, NMR is 
still capable to reproduce the correct dynamics of such 
pseudo-entangled states, since they behave like pure en- 
tangled ones [25] . 

Many works have been done about the dynamics 
of entanglement [24. 25 ,26 ,27]. It was always expected 
that the dynamics of the concurrence would give an 
asymptotic decay of the degree of entanglement, since 
it is the typical behavior of the decoherence processes. 
However, for some quantum states it was reported a 
curious phenomenon, named sudden death of entangle- 
ment, where the degree of entanglement vanishes at a 
finite time [27j . It has been recently measured in optical 
[25] and atomic [25] systems. 

Fig. 4 shows a NMR experiment of the sudden death 
of entanglement of the pseudo-entangled state = 
( 1 00) + |ll))/\/2, where its concurrence for such state 
is given by C 1 term, using the normalization procedure 
described in Sec. Ill on the relaxation of 23 Na. The 
continuous line is an exponential fit, included for com- 
parison (asymptotic decay). We can clearly see that the 
concurrence vanishes at a finite time just below 20 ms, 
which is the same order of the transverse relaxation 
time of the sample. 
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4.3 Discrete Wigner Function 

Another usual description to quantum states is the Wigner 
function formalism, which consists on the representa- 
tion of a quantum state in the phase-space. Wigner 
functions can be obtained from the so called "phase 
space point operators" [30] . They can be defined as the 
following expectation value: 



localized p = value, to become homogeneous (better 
seen in the inset, only for even values of momentum). 
This is an indication that there is no localization [34] , in 
the sense that it occurs in classically chaotic quantum 
maps. Finally, Fig. [5jd) exhibits the time-evolution of 
the interference term, W(7,p). One can clearly see the 
decoherence by comparing the initially red points (•) 
to the final blue ones (•). 



W(q,p) = Tr pA(q,p) 



(16) 



where the A(q,p) are the (Hermitian) phase space point 
operators and p is the density operator. Because we are 
working with discrete spin systems, it is more adequate 
the use of discrete Wigner functions |31j . Such a de- 
scription is made upon periodic boundary conditions for 
both, position and momentum (geometry of a torus on 
the phase-space). In what follows, we will use the den- 
sity operators (Eq[12|) relative to the relaxing pseudo- 
entangled state 2 -1 / 2 (| 00) + |H))- The corresponding 
phase-space point operators in matrix representation 
are shown on Appendix [Bj 

The decoherence of superimposed localized wave pack- 
ets is a well-known subject in the literature [32 , 33 . The 
Wigner function for the cat-state exhibits two well sep- 
arated positive peaks and a non-classical interference 
pattern between them, which oscillates between posi- 
tive and negative values. Such oscillations vanish due 
to the coupling of the system to the environment, lead- 
ing to a quantum-to-classical transition. Fig. [5] shows 
NMR results for the Wigner function of the cat-state 
( 1 00) + |11)) undergoing relaxation towards equi- 
librium. The duration of the experiment is about 70 
ms. Fig. 0(a) shows the discrete phase-space represen- 
tation of this state (left), as well as the equilibrium state 
(right). Notice that W(0,p) and W(6,p) are positive, 
whereas W(7,p) is an oscillating stripe between them. 
The oscillation of W(2,p) and W(A,p) are due to the 
use of periodic boundary conditions and the application 
of the discrete Fourier transform to obtain the momen- 
tum [31] . Such "images" are important for the calcula- 
tion of the marginal momentum from X^^o" 1 ^(?> P) t° 
guarantee it as a regular probability distribution. Fig. 
[5]Jb) shows the same plot as (a) for 3 and 12 ms, but 
now normalized according to the method described in 
Sec. III. We can clearly see the changes in the Wigner 
function caused by quadrupolar relaxation, in particu- 
lar the interference pattern. 

Time-evolution of the momenta distribution can be 
seen on the bottom of Fig. (He) . The initial distribution 
are the red bullets, (•), whereas the blue one (•) are the 
final values. Intermediate values are the open dots (o). 
From the plot we can see that momenta evolves from a 



5 Conclusion 

In this work a method for the normalization of the ini- 
tially pure part of a pseudo-pure NMR state undergo- 
ing relaxation was presented. The method allows fol- 
lowing directly the pure part of p, whereas keeping it 
a true density matrix, which is interesting for the the 
purpose of NMR quantum information processing ex- 
periments. We exemplified the method studying the re- 
laxation of initially pseudo-entangled states in which 
the phenomenon of sudden-death is observed, and the 
decoherence of Wigner functions in phase-space. The 
method is particularly useful for relaxation studies in 
the context of NMR quantum information processing, 
an area which still lacks experimental data. 
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A Redfield theory of pure quadrupolar 
relaxation spin 3/2 

The spin 3/2 system dynamics under the presence of relaxation 
processes can be described using the Redfield formalism for the 
deviation density matrix 1351 . The Hamiltonian that describes a 
I > 1/2 spin system in the laboratory frame, considering Zeeman 
and first order quadrupolar interaction can be expressed as 



(17) 



The first term describes the Zeeman interaction (with Larmor fre- 
quency of uil) and the second the static first order quadrupolar 
interaction (with quadrupolar frequency of u)q). Hence, for nu- 
clear spin 1 = 3/2 the eigenvectors of the Zeeman plus quadrupo- 
lar terms arc |3/2), 1 1/2) , | - 1/2) and | - 3/2), which can be 
labeled as 1 00) , |01), 1 10} and |11), corresponding to a two-qubit 
system. Here to simplify the expressions of the solutions (<r) of the 
Redfield theory found in 1181 we utilize A = a — cr^ . The charac- 
terization of the relaxation process is performed by means of re- 
duced spectral densities, which contain the information about the 
local field fluctuations |19II18| . Furthermore, when \ujq\ -C \^l\, 
the relaxation can be described by three reduced densities at the 
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Larmor frequency. Below we show A matrix elements in time 
function, 



B Matrix representation to phase space point 
operators 



A 12 (t) = A 12 e 



-t/r 12 



A 03 (t) = ^o 3 e" t/T12 , 



A)l (*) 



^23 (t) 



^02 (t) 



Al3 (t) 



(1+ 


e -t/T») 




2 


(1- 






2 


(1+ 


e -t/n) 




2 


(1- 


e -t/n) 



^01 



^02 ■ 



(1 



-t/r%\ 





2 


(1 


f e -*/u) 




2 


(1 


-e-'M) 




2 


(1 





^23 



^23 



-*/tqi 



-t/ T oi 



(18) 



(19) 



(20) 



(21) 



The phase space point operators are a set of Hermitian operators, 
which gives all marginal distributions obtained from the Wigner 
function 1311 . They are defined on a phase space grid of 2N X 2N 
points 



A{q,p) = —U q RV- p exp{iirqp/N), 



(28) 



where U is a cyclic shift operator in the computational basis 
(U\q) = \q + 1}), V is a cyclic shift operator in the basis related 
to the computational basis via the discrete Fourier transform and 
R is the reflection operator (R\q) = \N — q)). The phase space 
point operator form a complete orthonormal basis of the space of 
operators and 



A(q + £ q N, p + Z P N) = A(g,p)(-1)$» i+^p+^p 1 



(29) 



f° r £g>S<j = 0, 1. For this reason, it is clear that the N 2 phase 
- t / T 02 j (22) space point operators corresponding to the first N X N subgrid 
of the determine the rest. 

NMR 3/2 spin system have four eigenstates (computational 
basis, |3/2), 1 1/2) , | - 1/2) and | - 3/2)). Thus, the Hilbert space 
dimension is N = 4. The operators A(q,p) in matrix representa- 
t/TQ2 j (23) tion are shown below: 



^00 (t) 



-t/Ta +e -t/T I i 



-Ax> 



- e -2t/T 12+e -t/T^ 



(24) 



An (t) 



-2t/Ti 2+e -*/T 2 N 



-Au- 



A(0,p) 



A(1, P ) 



(25) A(2,p) 



(\ \ 

iP 
(-l) p 
\0 (-i)P / 



Tp/iV 



i2irp/N 



( 1 ^ 
?0 
(-l)P 

V (-t)P j 

I 1 \ 
iP 
(-l) p 
V {-i)P ) 



(30) 



(31) 



(32) 



^22 (*) 



o "Ax) + ^ o -^11" 



2 ^uu t 2 

( e -2t/T12 +e ,-VTl* 



(26) 



^33 (t) 



fe-t/TS+e-t/Ti'j ( e -2t/x 12+e -t/x 2 \ 
2 tA °° + ~ 2 ^U" 

( e -2t/xi 2+e -t/xi) 

.x. V L A. 



(27) 



where the elements A^i from right side are initial values of them 
(A = erg — Coo)- In addition, the characteristic times r are defined 
as function of the inverse averaged reduced spectral densities [18] . 
The characteristic times were determined from previous work as: 
roi = 4.6, T 02 = 4.7, Ti2 = 11.1, ri2 = 20.8 and r 2 = 23.8 in 
miliscconds. 



A(3,p) 



„i37vp/N 



( 1^ 

i p 
(-I)p 

V(-«) p o o oj 



(33) 



The projection of the density matrix p in these matrices (Eq |161 
yields a grid with only real values which represents the density 
matrix in discrete phase space. 
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Fig. 1 Calculated population values of the initial computational 
basis state |11) evolving under pure quadrupolar relaxation. The 
top and bottom plots represent two different types of time-fixed 
normalization to deviation density matrix, from which we ob- 
tain the pure part of the pseudo-pure state (see Appendix A for 
details). On the top scheme the matrix population become neg- 
ative after fO ms. On the bottom one populations are positive, 
but their initial values arc incorrect. 
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Fig. 2 Time evolution of the polarization and respective normal- 
ized populations for spin 3/2 nuclei in the initial state |11)(11|. 
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Fig. 3 Polarization for computational basis states encoded in the 
spin 3/2 nuclei. Also the initial and equilibrium limits to polar- 
ization arc shown. Notice that initial polarizations are e, whereas 
final ones are 21 e. Experimental data taken from previous exper- 
iments 1181 . 
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Fig. 4 Concurence decay of pseudo-entangled initial state 
2~ 1 ' 2 ( 1 00) + 1 1 1 ) ) _ The data were obtained from experimental 
deviation density matrices of the relaxation of 23 Na nuclei, nor- 
malized by the procedure developed in the section [3] The solid 
line is the exponential fitting, included for comparison. The de- 
polarization curve utilized to normalize the relaxing pseudo-cat 
state are the inside plot. 
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Fig. 5 The dccohcrence of the pseudo-entangled state ob- 
served from quadrupolar relaxation. In (a) and (b) are theo- 
retical and experimental discrete Wigner functions, where the 
bluc(red) squares represent the positive (negative) values and 
white squares are null values. Horizontal (vertical) axis is the posi- 
tion(momcntum) coordinates. The temporal evolution are shown 
in (c) for the discrete momentum distribution and in (d) for the 
oscillating term W(7,p). Red(Blue) dots are the initial(final) val- 
ues and the inside plot are the initial(rcd) and final(blue) momen- 
tum distribution for only even momentum. 



